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Abstract 

The Weibel instability is analyzed for quantum plasmas described 
by the Wigner-Maxwell model. For a suitable class of electromagnetic 
potentials, the Wigner-Maxwell system is linearized yielding a general 
dispersion relation for transverse electromagnetic waves. For a double 
Gaussian equilibrium with temperature anisotropy, the derived dis- 
persion relation generalizes the classical Weibel instability equation. 
More detailed analytical results are obtained for the cases of extreme 
temperature anisotropy and for a three-dimensional water bag dis- 
tribution. In all cases, quantum effects tends to weaken or suppress 
the instability. Applications are discussed for dense astrophysical ob- 
jects like white dwarfs and neutron stars as well as for tunnel-ionized 
plasmas with controllable perpendicular plasma temperature. 

1 Introduction 

Quantum plasmas have attracted renewed attention in recent years due to the 
ongoing miniaturization of ultra small electronic devices and micro mechan- 
ical systems [Ij, to the relevance of quantum effects for dense laser-plasmas 
and micro plasmas [2] and for dense astrophysical objects [3]. Quantum phe- 
nomena are relevant for these systems for a variety of reasons, the most usual 
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being the de Broglie wavelength of the charge carriers (electrons, positrons, 
holes and so on) becoming comparable to the characteristic dimensions of the 
system. Quantum ion-acoustic waves [I] , a quantum magnetohydrodynamics 
model [5J, shear Alfven modes in ultra-cold quantum magneto plasmas [6], 
quantum corrections for the Zakharov system [7]- [9] and nonlinear solutions 
for quantum magneto plasmas [TD] have been constructed. New quantum 
modes have also been identified for ultra-cold dusty plasmas [TT]— [TBJ, where 
quantum effects can be used for plasma diagnostics. The most recent de- 
velopments take care of spin effects in non relativistic quantum plasmas [17] 
as well as the associated magnetohydrodynamics equations [18], with possi- 
ble important applications for solid state plasmas as well as in the vicinity 
of pulsars and magnetars. In addition, there are the analysis of the non- 
linear instability of polaritons [19], of the dynamics of dark solitons and 
vortices in quantum electron plasma [20], of nonlinear interactions between 
intense circularly polarized electromagnetic waves and electron plasma oscil- 
lations [21], of the behavior of quantum diodes [22J, of nonlinear quantum 
dust-acoustic waves [23], of quantum ion-acoustic double layers [21], of lin- 
ear and nonlinear ion-acoustic waves in unmagnetized electron-positron-ion 
quantum plasmas [23], the construction of classes of solutions for the quan- 
tum Zakharov-Kuznetsov equation [26], of electron-acoustic solitary waves 
in dense quantum electron-ion plasmas [27J and of linear and nonlinear dust 
ion acoustic waves in ultra cold quantum dusty plasmas [28J. Finally, there 
are new experimental studies [2S] of weakly degenerate quantum plasmas in 
a gaseous regime (i.e. non solid state plasmas). A recent review on quantum 
plasma models and their range of validity can be found in |30j . 

In some systems [31], the ultimate influence of quantum mechanics in 
plasmas is the stabilization of some classically unstable mode, for sufficiently 
strong quantum effects. However, in the intermediate regime where quantum 
effects are not too intense but are nevertheless not negligible, there are situ- 
ations where unexpected quantum instabilities can arise. Examples on this 
are the quantum two-stream and three-stream instabilities [32] [M], showing 
unstable modes of pure quantum nature and no classical counterpart. In 
addition, unlike classical plasmas, there is no Penrose functional determining 
the linear stability properties of quantum plasmas [36] . These considerations 
points to the subtle role played in plasmas by quantum diffraction effects like 
tunneling and wave-packet spreading. Therefore, it is a relevant subject, to 
work out well known classical instabilities now in the context of quantum 
plasma models. In this perspective, the present work considers Weibel's in- 
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stability [37J. Weibel instability arises from temperature anisotropy in the 
equilibrium distribution function and is one of the fundamental instabilities 
of plasma physics. In the more recent years, it has been the central concept 
in several instances, like in fast ignitor scenarios [38], for particle accelera- 
tion and magnetic field generation in astrophysical settings [39l HQl HI], for 
collective non-Abelian Weibel instabilities in melting color glass condensates 
[12], in covariant relativistic scenarios [131 [44] . in electron-positron relativis- 
ts shocks [15] , with kappa and generalized (r, q) distributions [16] and in 
laser heated plasmas [4T] . 

The aim of this contribution is to get detailed information about the in- 
fluence of quantum effects on Weibel's instability. For this purpose, it is 
used the (kinetic) Wigner-Maxwell model, which is the quantum counterpart 
of the Vlasov-Maxwell system. In this context, the Wigner function plays 
the same role as the classical distribution function. Due to its mathematical 
complexity, frequently the electromagnetic Wigner equation [5] has not been 
taken as the basic tool in electromagnetic quantum plasmas. Rather, most 
works relies on the quantum hydrodynamic model [5J [33]. Here, however, 
the kinetic description is used in order to provide a easier comparison with 
previous results on Weibel's instability, which were constructed in terms of 
the Vlasov-Maxwell system. Nevertheless, notice that there are instances 
[4"8] where fluid models were also applied to Weibel's instability. Since the 
treatment is restricted to small amplitude waves, the electromagnetic Wigner 
equation can still provide meaningful results without too complicated ana- 
lytical difficulties. 

Sometimes, the Weibel instability is treated in conjunction with counter- 
streaming beams and/or ambient magnetic fields. However, as pointed out in 
|49j . in these cases it is more appropriated to talk about filamentation insta- 
bility. Indeed, the Weibel instability is prompted only by a single anisotropic 
system. So, here it is followed the original approach by Weibel [37], focusing 
only on the consequences of temperature anisotropy, but now allowing also 
for quantum effects. In addition, the analysis is restricted to non-relativistic 
systems. Notice that recently quantum effects were addressed for the fila- 
mentation (but not Weibel) instability [5U]. In this case, quantum effects 
have been shown to reduce both the unstable wave-vector domain and the 
maximum growth rate. 

This paper is organized as follows. In Section II, we write the disper- 
sion relation for small amplitude transverse electromagnetic waves in quan- 
tum plasmas, as derived from the Wigner-Maxwell system. Then it is as- 
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sumed a double Gaussian equilibrium distribution function where tempera- 
ture anisotropy is allowed. The resulting dispersion relation is analyzed for 
several limiting cases in Section III. Namely, there are considered the ultra- 
quantum case, the semiclassic case for small wave-lengths, and the semiclassic 
case for long wave-lengths In Section IV an equilibrium with extreme tem- 
perature anisotropy and a three-dimensional water bag model are considered, 
allowing for more detailed analytic results. In all cases, quantum effects are 
stabilizing. Section V discuss possible applications in the case of neutron 
stars and white dwarfs as well as for tunnel-ionized plasmas. As is shown 
in the following, the quantum effects are enhanced for larger density and 
larger temperature anisotropy. For neutron stars and white dwarfs, the high 
densities tend to enhance quantum effects. For tunnel-ionized plasmas, the 
densities are not so high, but the temperature anisotropy can be significantly 
enough. Section VI is reserved to the conclusions. 



2 Basic equations 

Consider a plasma composed of electrons (charge — e, mass m) and a neu- 
tralizing immobile ionic background. In terms of the Wigner distribution 
function / = f(r,v,t), the electron particle density n = n(r,t) and the 
current density J = J(r, t) are given by 



n 



J dvf , J = — e J dw f v . 



(1) 



All integrals are from minus to plus infinity unless otherwise stated. To 
proceed, it is necessary to work in terms of the electromagnetic potentials 
(0(r, t), A(r, £)), since the electromagnetic Wigner equation [SET] is written 
in terms of them and not the fields. In tangent space with coordinates (r, v), 
and time variable t, it reads 
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omitting the time-dependence of the several quantities and using summation 
convention in some terms. For the derivation of the electromagnetic Wig- 
ner equation in this form, it was assumed the Coulomb gauge, V • A = 0. 
Also, notice that the electromagnetic Wigner equation as in Eq. (3) of Ref. 
[5] was been written in phase space with canonical coordinates (q, p) and a 
time variable r = t. To derive (jSJ), it is necessary to transform using q = r, 
p = mv-eA and r — t, applying the chain rule, which implies, in particular, 
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In the formal classical limit, % — > 0, (j2J) reduces to the Vlasov equation, 
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(3) 
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(5) 



As a last remark, here we also have corrected a mistyping in one of the signals 
at the third term in the right-hand side of Eq. (3) of Ref. [5]. This mistyping 
produce no harm for the conclusions of this reference. 

Now take linear wave propagation along the Oz axis, with wave-vector 
k = kz, identically zero scalar potential and 



A = Aj_ exp(i[kz — uii\) 



(6) 



where Aj_ is a first-order object satisfying k • A^ = 0. This form for the 
vector potential is consistent with a transverse magnetic field. Also, the 
Wigner function expands as / = / + fx exp(i[kz — ut}), where / = /o( v ) is 
the equilibrium Wigner function, satisfying 



J dvf = n , J d\ f v = 



(7) 



where no is the ambient ion density and f\ is a first order perturbation. 
Accordingly, the current density is given by 



dv/iv 



(8) 
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omitting the exp(i[kz — ut]) dependence. 
In this context, Ampere's law yields 

{uj 2 - c 2 k 2 )A ± = — f dv/iv . (9) 

The homogeneous Maxwell equations are identically satisfied when working 
with the electromagnetic potentials, while Poisson equation is satisfied since 
it can be shown that there are no charge density fluctuations (/ dvfi = 0), 
as in the classical Weibel instability. Therefore, the only remaining equation 
is the linearized Wigner equation. Using (jSJ), the result is 

\ m ov J 

ev- Ai / hk hk \ 

= [fa{vx,Vy,v z -—)-f (v x ,Vy,v z + —)\, (10) 

for v = (v x ,v y ,v z ) and h = h/(2iv) being the scaled Planck's constant. 
Combining (j7j), (jUJ) and (fTUj) . assuming that the equilibria f are even func- 
tions of all velocity components also satisfying the condition fo(v x ,v y ,v z ) = 
fo(v v ,v x ,v z ) and doing an integration by parts, there follows the quantum 
dispersion relation for transverse waves (k • E = 0), 

2 2 2-2^ mUj2 p f, ( V l + V l \ v 

uj — uj — c k H ^ / av — - x 

p 2n h J \uj — fct^ / 

x ^/oK,%,^ + ^) - fo(v x ,v y ,v g - ^)J = 0, (11) 

where dj p = (n e 2 /(m£ )) 1 ^ 2 is the plasma frequency. 

Consider the specific case of the equilibrium Wigner function 

/»^er-(-^-?)-if). w> 

showing temperature anisotropy between parallel and perpendicular direc- 
tions. The temperatures T± and T\\ are measured in terms of energy units 
(Boltzmann's constant is Kb = 1). Then the dispersion relation ffTTj) develops 
into 

u?-M-<4(l + ^W Q )=0, (13) 
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where 

mv\\ ( r7 , uo hk . ry , uj hk .\ ^ A . 

w ^m\ z ^ + ^- z( Mr^)' () 

for v\\ = (2T\\/m) 1 ^ 2 and where Z is the plasma dispersion function. The func- 
tion Wq can be appropriately named a "quantum Weibel function". Apart 
from a scale factor, it is a centered finite difference version of the derivative 
of the plasma dispersion function. In the formal classical limit when H — > 0, 

(15) 

where the function W(u/(kv»)) is the same (Weibel) function as the one 
in Eq. (5) of Ref. [IE]. Hence, in this formal classical limit, the disper- 
sion relation ( 1131) reduces to the classical one, in accordance with the corre- 
spondence principle. Notice that quantum effects appears only through the 
non-dimensional parameter 

»=— , (16) 
mv\\ 

basically depending only on the longitudinal quantities k and v\\. 




3 Limiting cases 

We consider separately the ultra-quantum case, the small quantum parameter 
and wave-number case and the small quantum parameter and large wave- 
number case. 



3.1 Ultra-quantum case (H ^> 1) 

It is interesting to check the behavior of the dispersion relation in the case of 
very large or very small quantum effects. In this regard, an useful alternative 
(exact) expression for Wq is 

1 f d&'t 2 
W ° = 2V*J ({ -£)>-flV4- <17) 

for £ = uj/{kv\\). Assuming large quantum effects, so that H 2 ^> |c<j/(A;t>||)| 2 , 
it follows from (pU) that W Q ~ -2/H 2 , so that (JT3J) imply 

2 2,2 2/ 2m 2 vlT ± \ 

- = ck + < V-^k • (18) 
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There can be instability (u 2 < 0), provided there is also sufficient tempera- 
ture anisotropy, 

(19) 




However, since the right hand side of the last inequality is an increasing 
function of H, one conclude that quantum effects play an stabilizing role. 
Another way to derive this conclusion comes from rewriting (Tl9|) in terms of 
the instability condition 

-i+a+^V). (20) 




For increasing quantum effects, the critical wave-number k c becomes smaller, 
eventually dropping to zero. Notice also that ( 1201 and H 2 ^> 1 are conflicting 
conditions. A detailed but cumbersome analysis shows that both conditions 
can be satisfied only for hu) p /(mc 2 ) ^> 1, violating the non-relativistic as- 
sumption of the present model. 



3.2 Semiclassic case at small wave-lengths (H <C 1 and 

i«i « i) 

Now retaining only the first-order quantum correction, one get 

W Q = -1- £Z(£) + £ (2 + 3£Z(£) - 2e - 2eZ(0) + 0(H 4 ) . (21) 

For such small quantum effects and also taking small wave-lengths so that 
|f | < 1 and Z{£) ~ iy/v, there follows from (EED that 

W Q ~-l-i^+^-, (22) 
while the dispersion relation (T3~3|) produces 

ikv\\T\\ {T ± , #\ c 2 k 2 \ , , 

In the derivation of (1231 . it was used \uj 2 \/(c 2 k 2 ) = \i 2 \v 2 /c 2 1, also 
consistent with the non-relativistic approximation. Assuming positive wave- 
numbers, it is apparent from ( 123]) that purely growing waves (frequencies 
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composed only by a positive imaginary part) can exist for sufficiently large 
temperature anisotropy, and that this necessary anisotropy becomes larger 
for increasing H. Specifically, for small quantum effects and wave-lengths, 
there will be a purely growing wave only if 

71 1 H 2 

showing the need of extra anisotropy for instability, due to quantum effects. 
In addition, in this combined regime of small quantum effects and wave- 
lengths, it can be proven that the unstable wave-numbers are restricted to 



C 2 ( ! + I^L 



k , 2 _ -pv-^-ii (25) 

Once again, the stabilizing nature of (now small) quantum effects is apparent, 
since the allowable unstable wave-numbers occurs for a smaller range for 
increasing H 2 . 



3.3 Semiclassic case at large wave-lengths (H <C 1 and 

Kl » i) 

From the expansion (|21|) and using Z(£) ~ — 1/£ — l/(2£ 3 ) when |£| ^> 1, for 
large wave-lengths, it follows the dispersion relation 

^-^-0,^1 + ^(1-^=0. (26) 

For |u;| ck, ( f26l) yields the purely growing mode 

Notice that the growth rate becomes smaller for larger H. Also, in view of 
( 127)1 . the condition |f | > I can be attained only for T ± > 7j|/(l - H 2 /A). 
Equations (|26j) and (1271) are formally the same as Eqs. (8) and (9) of Ref. 
|48j . making the replacement 71 — > 71(1 — 77 2 /4), indicating the need of 
extra temperature anisotropy, in view of quantum effects. 
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4 Toy models and more detailed analytical 
results 



Anisotropic Gaussian distributions are amenable to analytical results only 
for limiting situations. This Section consider some toy models which behave 
more friendly in this respect. As a first example where full analytical results 
are available, consider the equilibrium Wigner function 

/o= 2^^ )eXP ("2^ ( ^ + ^)' (28) 

which can be viewed as a distribution with extreme temperature anisotropy 
(Tj| — > 0). Inserting into (iTTj) and proceeding as before, it results 



2 \ p Am 2 



> 2 + c 2 k 2 - — — 2 + 



p 



4m 2 m 



• (29) 



One of the roots is unstable {u 2 < 0), provided 



showing, once again, stabilization due to increasing quantum effects. 

The classical transverse Weibel instability has sometimes considered in 
terms of water bag distributions [38], S3 E2J- As a second example, also 
amenable to detailed calculations, take the following three-dimensional water 
bag [4]5] equilibrium, 

fo = (0(v x + v ± ) - 6(v x - v ± )) {6(v y + v ± ) - 6(v y - v±)) x 

O t/j^1/jj 

x (e{v z + v {l ) -6{v z -v l{ )) , (31) 

where 9 is the Heaviside function and in this context v± and v\\ are related 
to dispersion of velocities in the perpendicular plane and along the Oz axis, 
as before. Then the dispersion relation (iTTj) yields 



22 muj 2 v\ I J 1 - {kv\\ - hk 2 /(2m)) 2 * 
6hk Vll 11 [u 2 - (k Vll + hk 2 /(2m))\ 



- l - < - ^ = 1- , I , : , Z 2 > (32) 
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not in polynomial form as in the formal classical limit for the corre- 
sponding equilibrium. However, some analytical results are still available. 
Assuming purely growing instabilities with a growth rate 7, so that 00 = 1^, 
and also disregarding 7 2 at the left-hand side of (132]) assuming 7 2 -C cu 2 , 
there follows 



7 = coth ] ^— 2 v - - k 2 v\ - — — . 33 

m \ m(jjpVj_ J " Am z 

Notice that (133]) is not valid in the formal classical limit fi = because then 
the associated growth rate would not be small. Rather, in this limit one has 
to Taylor expand the right-hand side of ( 1321) . and then the results from [19] 
are recovered. 

To proceed, it is convenient to adopt the following rescaling, 

7 1 kv \\ fj huJ P v ± v \\ r oA\ 

7 = — , k = — - , h = — 2 > v ± = — > n = — » ( 34 ) 

ujp ujp mvf, c c 

so that ( 1331) becomes 

7 2 = M 3 coth (^{k 2 + t; 2 )^ - ¥ - ^ . (35) 

It is possible to estimate the maximum wave-number for instability, using 
coth(£) — l/£ for |^| 1. Assuming that this expansion is valid, one will 
conclude from (|35|) that 7 2 > provided 

k 2 <k 2 m ^^- 2 (i(H 2 vl - 4) 2 + 1^!1]V2 _ ^2,2 _ 4 ^J (36) 

From (|36|) . it can be shown that there will exist some unstable mode (k^ > 0) 
if and only if there is sufficient temperature anisotropy, v± > \/3v\\ , indepen- 
dently of the strength of the quantum effects, and in accordance with the 
scaling found for classical plasma [15] . However, it can be deduced from 
that the critical wave-number shrinks to zero as H —>■ 00. 



5 Applications 

Consider now the case of anisotropic Maxwellian equilibria. For semiclassic 
and small wave-lengths conditions (Section III. 2), one can use (1231) to show 
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that the maximal value of the growth rate happens at k = k m = k c / v3, with 
k c given by (J2SJ) • The corresponding maximal growth rate is then 

( § \ 1/2 / 7|| V /2 ^ll gVZj - Iff! ^ 
V277T/ y \mc 2 ) T\ /_ T,h 2 Lol\ ' 



1 + 



12T,fmc 2 



The calculations leading to (1371) were semiclassic in the sense that the pa- 
rameter H = hk/mv\\ is taken as small. However, there can be significant 
deviations from the classical expression, coming from the h term in the de- 
nominator of (|37p . provided there is sufficient temperature anisotropy. These 
deviations are also enhanced by large densities. For instance, consider white 
dwarfs and neutron stars, with typical values n ~ 10 32 m~ 3 , T± = 10 7 K, 
T\\ = Tj_/100. The origin of the temperature anisotropy can be, for in- 
stance, the propagation of a shock wave. For these parameters, for hydrogen 
plasma, one finds that the quantum corrected maximal growth rate is about 
11% smaller than the classical one. However, these calculations have to be 
taken with care, because they suppose that H 2 from (fT5T) at k = k m and 
|£| — 1m/(k m v\\) are small quantities. For the chosen parameters one get 
H 2 ~ 0.41 and |£| ~ 0.37. For larger densities, H 2 would be even greater. 

Another interesting system where quantum corrections for Weibel insta- 
bility can be significant are tunnel-ionized plasmas with negligible longitudi- 
nal temperature and where the perpendicular temperature can be controlled 
by a varying laser polarization. It has been argued [S3] that the Weibel in- 
stability could be a mechanism for further increase of Tji with time. For a 
typical value [53J of T\\ ~ leV, one find 

T±h 2 u 2 „. 

=—5- = 2.3 x 10" 34 n T ± , 38 

UTfmc 2 ' K ' 

using SI units for n and T±. Although it is not easy to get large values 
for the quantity at the right-hand side of (138|) . the fast progress in next 
generation intense laser-solid density plasma interaction experiments can be 
such that quantum effects stops the Weibel instability, specially for ultra- 
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dense systems. For the largest densities [29], [Mj feasible now, no ~ 10 
and for T± ~ lOOeV, the right-hand side of (138]) has already a significant 
value of the order of 25. Another way to get even larger quantum effects is 
the use of smaller values of T». The results of this section, however, have to 
be taken with care, since they deal with dense plasmas which should be more 
properly treated by Fermi-Dirac statistics. 
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6 Conclusions 



In general terms, quantum effects produces smaller growth rates and smaller 
ranges for unstable wave-numbers, in the case of equilibria with distribution 
functions anisotropic in temperature. Some applications in astrophysical 
scenarios and in tunnel-ionized plasma were discussed. Eventually, in the 
ultra-quantum case, the unstable region shrinks to zero. We can understand 
this result in an heuristic way as follows. Due to wave-particle spreading 
and tunneling, quantum effects tends to enhance the dispersion of parti- 
cles in phase space. This corresponds to an effectively smaller temperature 
anisotropy, or thermalization, so that the original ratio T±/T\\ have to be 
greater to produce the same instability results as in classical plasma. Sim- 
ilar spreading in phase space also occurs in the case of quantum corrected 
Bernstein- Greene-Kruskal modes [55] - However, as mentioned in the Intro- 
duction, sometimes quantum effects can also gives unexpected enhancement 
of plasma instabilities. Therefore, it is interesting to pursue this trend, look- 
ing at the behavior of additional well-known classical plasma instabilities, 
in the context of quantum plasma models. Also, an important issue is the 
inclusion of relativistic and spin effects. 
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